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The energy of fluctuating electromagnetic field is investigated for the thermal Casimir force acting between parallel plates 
made of real metal. It is proved that for nondissipative media with temperature independent dielectric permittivity the energy 
at nonzero temperature comprises of the (renormalized) energies of the zero-point and thermal photons. In this manner 
photons can be considered as collective elementary excitations of the matter of plates and electromagnetic field. If the dielectric 
permittivity depends on temperature the energy contains additional terms proportional to the derivatives of e with respect to 
temperature, and the quasiparticle interpretation of the fluctuating field is not possible. The correlation between energy and free 
energy is considered. Previous calculations of the Casimir energy in the framework of the Lifshitz formula at zero temperature 
and optical tabulated data supplemented by the Drude model at room temperature are analysed. It is demonstrated that this 
quantity is not a good approximation either for the free energy or the energy. A physical interpretation of this hybrid quantity 
CNj ' is suggested. The contradictory results in the recent literature on whether the zero-frequency term of the Lifshitz formula for 
the perpendicular polarized modes has any effective contribution to the physical quantities are discussed. Four main approaches 
to the resolution of this problem are specified. The precise expressions for entropy of the fluctuating field between plates made 
of real metal are obtained, which helps to decide between the different approaches. The conclusion is that the Lifshitz formula 
y—( ' supplemented by the plasma model and the surface impedance approach are best suited to describe the thermal Casimir force 
CO between real metals. 

PACS numbers: 12.20.Dc, 42.50.Lc, 65.50.+m 
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The Casimir effect is a rare macroscopic manifestation of the zero-point electromagnetic energy. It results from 
the alteration of the zero-point spectrum by the material boundaries (see the original Casimir's paper jjj and ex- 
tensive reviews ) . I n recent years considerable attention was paid to the precision measurements of the 
Casimir force between metallic surfaces fip|,p| JTc|Jll|Jl^ , p^ . p^^5| . The results of these measurements were used in 
Refs. |l6||L^jn||l^^0[ to constrain the hypothetical interactions predicted by many extensions to the Standard Model 
and also in nanotcchnology . This called for new theoretical investigation of the Casimir force with allowance made 
for the realistic boundary properties, i.e., surface roughness, finite conductivity of a metal and nonzero temperature 
(see Ref. || for the review). Also the combined effect of these factors has attracted considerable attention. 

In this paper we consider the correlation between energy and free energy for the Casimir force acting between two 
plane parallel plates made of real metal of finite conductivity kept at nonzero temperature. At zero temperature the 
, influence of the finite conductivity of the boundary metal onto the Casimir force was examined in Refs. |2^ , ^3| , ^ , p5| 
on the basis of the Lifshitz theory using the optical tabulated data supplemented by the Drude model. In this 
approach the Casimir energy was represented by an integral with respect to continuous frequency, as it is at zero 
temperature. However, the optical tabulated data and the value of the relaxation parameter of the Drude model at 
T = 300 K were substituted, which are actually temperature dependent. By this means, in fact, some hybrid quantity 
was computed, which is different from both energy and free energy, with no clear relevance to either of them. When 
it is considered that the Casimir energy between two plates is generally used to calculate the Casimir force in the 
experimental configuration of a sphere (spherical lens) above a plate, the resolution of this issue is of great interest. 
Below is shown the relationship between energy, free energy and a quantity computed in Refs. |2^ , p3] , p4| , p5| . 

Investigation of the correlation between energy and free energy for the thermal Casimir force has also assumed great 
importance in connection with the contradictory results obtained by different authors when applying the Lifshitz 
theory to real metals H,|^,||j2|j3^,|l|,||,|||3|,|| . The contradictions arise on whether or not the zero-frequency 
term of the Lifshitz formula for the perpendicular polarized modes of electromagnetic field contribute to physical 
quantities. There are four main approaches to the resolution of this issue in the recent literature. 

a) According to Refs. [ p6||35| ] based on the immediate application of the unmodified Lifshitz formula, the zero- 
frequency term of this formula for the perpendicular polarized modes is equal to zero in the case of real metals 
described by the Drude model (remind the reader that for ideal metals the reflection coefficients for both polarizations 
are equal to unity at zero frequency, and hence the zero-frequency term for both modes is not equal to zero). This 
approach leads to the conclusion that thermal corrections to the Casimir force are large, negative, and linear in 
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temperature at small separations, and the asymptotic Casimir force between real metals at large separations is two 
times smaller than for the case of an ideal metal (with no regard for the particular value of the conductivity) 133] . 



b) From the standpoint of Refs. 28 32], the perpendicular polarized modes give a nonzero contribution at zero 
frequency. To find it, a special modification of the zero- frequency term of the Lifshitz formula was proposed (found 
by analogy with the prescription of Ref. ]36| ] for an ideal metal but not coinciding with it). 

c) According to Refs. |3^,^l|, the perpendicular polarized modes also give a nonzero contribution at zero frequency. 
For real metals it is, however, the same as for ideal metal, i.e., the reflection coefficients for both modes arc equal to 
unity at zero frequency. What this means is for real metals the same modification of the zero-frequency term of the 
Lifshitz formula is made as for ideal metals This approach leads to linear (although positive) thermal corrections 
to the Casimir force at small separations and to the absence of any finite conductivity corrections for real metals 
starting from moderate separations of several micrometers regardless of metal quality ]32j . 

d) Finally, according to the approach of Refs. p7| , p8|j29| both modes with the parallel and perpendicular polarizations 
do contribute to the zero-frequency term and this contribution can be calculated by the substitution of the plasma 
model dielectric function into the unmodified Lifshitz formula. The same conclusion is obtained in Ref. [fj3| on the 
basis of the surface impedance approach. 

Thus, at the moment there is no agreement in the theoretical literature as to the description of the thermal Casimir 
force between real metals. To gain a more complete understanding of the present state of affairs, in Ref. ||J] the 
thermodynamical argument was exploited. According to Ref. J34|, the approaches a) and c) do not conform to the 
requirements of thermodynamics as they lead to the negative values of entropy and violation of the Nernst heat 
theorem. Although the qualitative conclusions of (3^] are quite correct, the quantitative calculations are incomplete 
as they do not take into account the entropy of real photons. The precise expressions for the energy and free energy 
of the fluctuating electromagnetic field found below are used to obtain the quantitative behavior of entropy as a 
function of surface separation distance and temperature. The obtained results confirm the conclusion of Ref. ]3^] that 
the approaches a) and c) are not compatible with thermodynamics. They also give the possibility to compare the 
approaches b) and d) in order to decide between them. 

The paper is organized as follows. In Sec. II the main notations are introduced and the case of nondissipative 
condensed media separated by a gap is considered, with the media described by a temperature independent dielectric 
permittivity. It is proved that in this situation one can introduce photons as quasiparticles due to the collective 
elementary excitations of condensed matter and electromagnetic field. As a consequence, the energy at temperature 
T defined via the derivative of the free energy with respect to T comprises of the (renormalized) energies of the 
zero-point and thermal photons. In Sec. Ill the energy and free energy of the fluctuating electromagnetic field are 
considered on the basis of the Lifshitz theory and the plasma model. In Sec. IV the applicability of the Drude model 
for the calculation of the thermal Casimir force is discussed. The energy and free energy are found in the case of a 
metal described by the Drude model. The alternative approaches to this problem available in literature are analysed 
and compared. The physical sense of the energy at temperature T appears to be more complicated than in the case 
of the plasma model. It is shown that energy contains additional terms depending on the derivatives of the dielectric 
permittivity with respect to temperature. In Sec. V the entropy for the thermal Casimir force acting between real 
metals is calculated precisely for both plasma and Drude dielectric functions. Sec. VI contains conclusions and 
discussion. 

II. PHOTONS BETWEEN PLATES AS ELEMENTARY EXCITATIONS 

We consider the configuration of two semispaces (thick plates) with frequency-dependent dielectric permittivity 
s(uj) restricted by parallel planes and separated by an empty space with distance a between them at a temperature 
T. This is a system in thermal equilibrium. The free energy per unit area is given by the well known Lifshitz formula 
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Here A|i t ±(u), k±) are the quantities having zero values on the photon eigenfrequencies permitted between the plates 
by the boundary conditions (indices ||, _L label two independent polarizations, and k± is the modulus of a wave vector 
in the plane of plates) 



A||,x(^ n) fc±)=0. (2) 
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They can be expressed in terms of reflection coefficients on the imaginary frequency axis 

A||,xte,*:±) = l-rJ i± (6,Aj.)e- a ^, 



where 



with the notations 
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In Eq.(|l|) k B is the Boltzmann constant and = 2nlkBT/h, where Z = 0, ±1, ±2, . . . , are the Matsubara frequencies. 
As seen from Eq.(^), the quantities Aim are normalized in such a way that the free energy (Q) tends to zero for the 
infinitely remote plates (a — > oo). The details of the renormalization procedure can be found in Refs. P,p5|. 

In this section we consider nondissipative media, which is to say that e(oj) is a real function. It will be assumed also 
that at a given frequency e does not depend on temperature. Both conditions are satisfied, e.g., for metals described by 
the plasma model or for dielectrics with a constant dielectric permittivity (the case of dissipative media is considered 
in Sees. IV, V). Under these conditions we find the simple quasiparticle interpretation for photons between plates and 
the expression for the energy of the fluctuating field at a temperature T. 

According to thermodynamics, energy at arbitrary temperature is given by 



E{a,T) = -T 



a 8 F E (a,T) 
dT T 
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where the free energy is defined in Eq. ([!]). 

Taking into account that the term of Eq.(0) with I = is linear in temperature and the quantities Aim are even 
functions of I one obtains 



E(a,T) = —?—Y^j Q k ±dk ± — [lnA||te,^) +lnA ± te,^)] . 
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Let us next use that Aj| j_ depend on temperature through the Matsubara frequencies only, so that d/dT 
(£l/T)d/d£i. Thus, energy per unit area is given by 
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E(a,T) = — / fcj.dfcj.^-[lnAyte > fc J .)+]nA ± te J fcj.)]. 
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We consider now the interpretation of energy at temperature T in terms of elementary excitations. In the case of 
the nondissipative media under consideration the photon eigenfrequencies are real and the nonrenormalized energy 
of equilibrium fluctuating electromagnetic field in the system comprises of the energy of zero-point fluctuations and 
Planck's protons |39] 



Identically, Eq. ([)]) can be rearranged to give 
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This quantity is evidently infinite. The renormalized value of the sum over the eigenfrequencies wf'^ can be 
calculated by the use of the argument theorem |flp|j6| j2qj40|| . 
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where the integration path C in the plane of complex to is shown in Fig. 1, and the normalized quantities A|| j_, 
having zero values on eigenfrequencies, were substituted defined by Eqs.(|])-(||), with changed by —iuj. Note that 
the function uj coth[huj / (2ftBT)] has poles at the imaginary frequencies uji = I = ±1, ±2, . . . , where £/ are the 
Matsubara frequencies (it is, however, regular at luq = 0). Because of this, the integration along the imaginary axis 
involves semicircles about these poles. Integration along a semicircle, whose radius extends to infinity, makes zero 
contribution to the right-hand side of Eq.dnl). As a result Eq.(nil) leads to 
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Here prime denotes the derivative with respect to to and we took into account that Aim are even functions of lu, 
so that their derivatives are odd ones. This property leads also to the zero value of the seemingly pure imaginary 
integral in the right-hand side of Eq. (|T^). 

Substituting the right-hand side of Eq. (|l^ ) into Eq .([ic| ) instead of a nonrenormalized sum, we obtain the renormal- 
ized energy at a temperature T coinciding with Eq.(|8|) derived from the thermodynamical definition (||). In such a 
manner we have proved that at certain conditions the thermodynamical energy at equilibrium is given by the additive 
sum of the contributions from the zero-point fluctuations and Planck's photons. The renormalization of both quanti- 
ties reduces to the subtraction of the contribution of a free space with no plates. As a consequence, in the absence of 
dissipation, photons between plates can be considered as some kind of quasiparticle excitations in the system of the 
electromagnetic field interacting with the matter of plates. The simple example of this situation is given by metals 
described by the plasma model. 



III. ENERGY AND FREE ENERGY FOR THE THERMAL CASIMIR FORCE IN THE FRAMEWORK 

OF THE PLASMA MODEL 

The considerations of the previous section can be illustrated by the dielectric function of the plasma model 
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e ( w ) = l-% £ (^) = 1 + ^, (13) 

where lu p is the plasma frequency. This dielectric function is real and its parameter does not depend on temperature. 
The use of the plasma model to calculate the thermal Casimir force corresponds to the approach d) described in 
the introduction. The free electron plasma model works well for frequencies of visible light and infrared optics. 
It is common knowledge that the dominant contribution to the Casimir effect comes from the range around the 
characteristic frequency lu c = c/(2a). Thus the plasma model is applicable in the a-range from a few tens of nanometers 
to around a hundred micrometers. 

For the sake of convenience, we introduce the dimcnsionless variables 

f=M y 2 =4a 2 (k 1 _ 2 + ^) (14) 



c 



in terms of which the plasma dielectric function takes the form 



~ 2au) v 



e (^) = l + 1 §, u p = ^. (15) 



In terms of the variables £, y the Lifshitz formula (hj) can be rewritten as 
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where 



(17) 



and the reflection coefficients are 

rf\(ihv) = \ 
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By virtue of the fact that e depends on temperature through the Matsubara frequencies only, the zero-frequency 
term of Eq. ( [Tif ) (I = 0) does not contribute into the energy (0) [compare with Eq. J7l)]. It is notable also that in the 
special case of the plasma model the perpendicular reflection coefficient from Eq. (jig) is given by 



rl(£i,y)=rl(y) = 



y 



UJ 2 



(19) 
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i.e. it is frequency- and temperature-independent for any I. By this reason its derivative with respect to temperature 
does not contribute to energy @. As a result, in the framework of the plasma model the energy per unit area at a 
temperature T, calculated by Eqs. (ra), (|l6|), takes the form 



1 rp oo 
1=1 



if [lnA||(6,6) + lnA ± (6,e0 



Ana 2 . 
1=1 J ^ 



ydy- 
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(20) 



This equation is convenient for numerical calculations. 

Let us now compare the values of energy at temperature T given by Eq. 
values of energy at zero temperature given by [p|rj, |37| , |38|] 



]) and free energy of Eq. (jlj) with the 



E(a,0) 



he 



/ ydy 



lnA||(6y)+]nAj.(£,iO 



327r 2 a 3 

The calculational results at T = 300 K for the case of Al with |tl| 

lo p = 11.5 eV= 1.75 x 10 16 rad/s 
are shown in Fig. 2. In this figure the dimensionless ratios 

EP l {a,T) Fg(a,T) £ p '(a,0) 



\EP l (a,0)\' |J5P'(o,0)|' \EP l (a,0)\ 



-1 



(21) 



(22) 



(23) 



are plotted by the solid lines 1,2 and dashed line, respectively, as the functions of the surface separation. The energy 
at zero temperature E pl (a,0) is computed by Eq. ( ^l| ) where the plasma dielectric function given by Eq. ( |l3| ) is 
substituted. It is clearly seen, that at smallest separations all three quantities (energy at T = 0, energy and free 
energy at T = 300 K) have approximately equal values. With an increase of the separation distance the modulus of the 
relative energy at temperature T decreases to zero limiting value while the modulus of the relative free energy increases. 
Note that the limiting cases of small and large separations can be simultaneously considered as the limits of low and 
high temperatures, respectively, if one compares with the so called effective temperature ksT e ff = Tiuj c = hc/(2a) 

iis 

The asymptotic behavior of energy and free energy at small and large separations (low and high temperatures) 
in the case of the plasma model can also be investigated analytically. As was proved in Ref. J2S[| , one can expand 
Eq. (16) in powers of a small parameter A p /27ra, where X p is the plasma wavelength, and in a contribution, depending 
on temperature, it would suffice to preserve the first power only. The result valid for all a > X p is 
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(24) 
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where t = T e ff/T. The quantity E pl (a,0) is the energy at zero temperature. Its expansion in powers of Xp/lim 
can be found in Refs. (here the result up to fourth order should be used in order to get sufficient accuracy at 

smallest separations). 

From Eq. (^J) the required asymptotics follow. At small separations (T <C T e ff) one obtains 



na 



T 



T, 



eff 



45C(3) 



1 + 2^ 



T 



T, 



eff 



(25) 



where £(z) is the Riemann zeta function. 

Applying the thermodynamical definition (^|) to Eq. (p5|), we obtain the low-temperature asymptotic of energy 



E pl (a,T) = E pl (a,0) + ^f^ 

OTTCL 



2E 
Tin 
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30C(3) 
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T, 
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(26) 



This asymptotic expression is obtained also from Eqs. (|J) or ( p0| ) by the use of the Abel-Plana formula (see Ref. p£ 
where similar calculations were performed). 

In the opposite case of large separations (T ^> T e ff), Eq. ( ]24| ) leads to the main contribution of the form 



*g(a,T) = -^C(3) fl-^ 



(27) 



By virtue of Eq. (^) the asymptotic value of energy is E pl (a, T) = 0. If one wished to have a more exact asymptotic of 
energy, the next (exponentially small in T/T e //) terms omitted in Eq. (27) should be taken into account or Eq. ( |20| ) 
should be used. In both cases the result is one and the same 
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Comparison of the numerical calculations presented in Fig. 2 with calculations by the asymptotic formulas of 
Eqs. (p5|)-(p8|) shows that the asymptotic of small separations works well within the separation range A„ < a < 
2 — 3/xm, and the asymptotic of large separations is applicable for a > 5 fim. In the transition range, Eqs. (pi]), j20| ) 
should be used to calculate the values of the free energy and energy for the thermal Casimir force in the framework 
of the plasma model. 

If we consider the limit up —> oo (Xp —> 0) in Eqs. ( 24 )— (|2q), the results for ideal metal are obtained. 



IV. DIFFERENT APPROACHES TO THE CALCULATION OF ENERGY AND FREE ENERGY IN THE 

FRAMEWORK OF THE DRUDE MODEL 

Let us now consider metals described by the Drude dielectric function 
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sH = l- , e(l0 = i + _^, (29) 



where 7 is the relaxation parameter. In terms of a dimensionless frequency introduced in Eq. (14), the Drude dielectric 
function along the imaginary axis is 

As was noticed in the introduction, there is no agreement in the recent literature regarding the use of the Drude 
model in the framework of the Lifshitz theory. Because of this, it is appropriate to reexamine the applicability of 
the Drude model in the context of the thermal Casimir force. The Drude model, as opposed to the plasma model, 
takes into account the phenomenon of volume relaxation. In reality, this phenomenon plays a role in the domain 
of the normal skin effect where the mean free path I of the electron is much less than the penetration depth of the 
electromagnetic oscillations into a metal S and the mean distance v/uj traveled by an electron in a time l/2w of the 
period of the electromagnetic field p3|,E4| . For most of metals at T = 300 K the domain of the normal skin effect 
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extends from the quasistatic fields (where e is pure imaginary) to the frequencies of order 10 12 rad/s. What this means 
is that the Drude dielectric function has a direct relationship only to plate separations 0.1cm < a < 1km such that 
the characteristic frequency lu c — c/(2a) belongs to this domain. However, at so large separations the Casimir force 
is extremely small and is of academic interest only. 

For higher frequencies, depending on which metal is considered, the anomalous skin effect (8 <C I, 8 -C v/lu) or 
relaxation region (v/lu <C I <C 8) occur. Here the volume relaxation described by the parameter 7 is not significant, 
but, in general, the space dispersion gives an important contribution. Note that in the domain of the anomalous 
skin effect (it extends up to around 7 x 10 13 rad/s) metal cannot be described by either the Drude model, given by 
Eqs. (p9|), (|30"|), or by any dielectric function depending only on frequency as well. 

On further increase of frequency, the transition to the infrared optics occur, where v/lu <C 8 <C I (or to the 
"extremely anomalous skin effect" if we use Casimir's terminology ]43|]). In this domain the volume relaxation does 
not play any role. In the semiclassical theory of AC conductivity e is practically real, signifying no dissipation of 
the electromagnetic energy within the metal jis)] . Because of this, the plasma model is realistic if the characteristic 
frequency lu c belongs to the domain of the infrared optics (see Sec. III). This domain extends to the frequencies of 
around 2 x 10 16 rad/s and for higher frequencies is followed by the domain of the ultraviolet transparency of metals. 
However, some interelectron collisions and a scattering on the surface lead to a small imaginary part of £ in the domain 
of infrared optics Eq] as is demonstrated by the optical tabulated data for complex refraction index |4l|] . These data 



are often used to find the values of e(i(,) along the imaginary axis through the dispersion relation f]|, l22| , [23 24 25f| . By 
way of example, for Al the optical data for u > 6.08 x 10 13 rad/s are tabulated ftjj . 

At the same time, the existence of the anomalous skin effect domain, where the concept of s(lu) is not applicable, 
is usually ignored, and the optical tabulated data are theoretically extended into the domain of lower frequencies by 
means of the Drude dielectric function 41 1. This is needed to compute the dispersion integral from zero to infinity. 
The values of s(i(.) obtained in such a manner by means of the dispersion relation and extended tabulated data are 
satisfactory up to £ ~ 10 15 rad/s with e(i£) obtained by the immediate substitution of the imaginary frequency into 
the Drude model according to Eq. ( p^ ) with no use of the dispersion relation. This suggests that the Drude model can 
be applied for the calculation of the Casimir force within a micrometer domain a > 0.4 fim in parallel with the plasma 
model. It should be particularly emphasized, however, that the application of the Drude model in the domain of 
infrared optics is physically unjustified as the volume relaxation is absent in this domain (below we call into question 
also the possibility to substitute the Drude dielectric function into the zero- frequency term of the Lifshitz formula). 

In contrast to the case considered in Sec. II, the Drude metals are dissipative media, described by the complex 
s(lu). At a given frequency e depends explicitly on temperature through the relaxation parameter 7. Because of this, 
the energy of the equilibrium fluctuating electromagnetic field cannot be presented any more in the simple form of 
Eq. (JsJ) . In accordance with the thermodynamic equality , additional terms appear in the right-hand side of Eq. (|^) 
containing the derivatives duj^ n /dT . 

Substituting Drude dielectric function ( 30|) into the Lifshitz formula for the free energy ( |l6|) and using the definition 
(^) of energy at temperature T, one obtains 



E D (a,T) = - 



k B T 2 8 
167ra 2 df 
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9r\\(ji,y) T dr\\(ji,y) 



dj dT 



x dr±(ii,y) _9r x (6,y) ^7 
6 - +T 3 7 df 



Here the zero-frequency term of Eq. ( |i6| ) is separated because there is disagreement in recent literature on whether 
or not it contributes to the Casimir energy and force. The immediate consequence of Eqs. (||), (|l^), (^) [approach 
a) described in Introduction] is pq] 



dyyln (l 



-C(3). 



This result is given by the parallel modes only, while the perpendicular modes do not contribute. 

The special modification of the zero-frequency term of Eq. (|l6|) proposed in p3| [approach b)] leads to 



/>oo 

/( b) (a,T) = -C(3)+ / dyy\n[l-rl(y,y)e-y] 
Jo 



(32) 



(33) 
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The two contributions in the right-hand side of Eq. (|33|) are given by the parallel (perpendicular) modes, respectively. 
If for real metals the same prescription is used as for ideal metal [approach c)] , one obtains ]30|,pl[ 



dyy]n(l-e- y )=-2{(3), 



(34) 



where both polarizations lead to equal nonzero contribution. Evidently, in the framework of the approaches a) and 
c) the zero- frequency term is temperature independent and does not contribute to the energy ( |3l"l ) . In the framework 
of the approach b) there is only a fair contribution due to the dependence of r±(y, y) on 7(F) in Eq. ( |33|). 

Note that Eq. (31) is in direct analogy to Eq. (f20|). The additional terms which are present in Eq. (|3l|) take into 
account the explicit dependence of the dielectric permittivity on temperature through the relaxation parameter. This 
equation is convenient for the numerical calculations. 

Before performing the calculations, let us give the approximate expressions for both free energy and energy which 
allow one to compare the results obtained in the framework of the Drude and plasma models. For this purpose we 
expand Eq. ( |l6| ) in powers of a small parameter 7 /u> p preserving the first-order term only (for Al at T — 300 K 
7 = 0.05 cV = 7.6 x 10 13 rad/s, so that for lower T 7/wp < 0.004). The coefficient near this term can be computed in 
the zeroth order in a small parameter a — Xp/Ana — l/u> p . The result is 



FE(a,T)=Fg(a,T) 



7 k B T 



POO >t 

/ (a ' b ' c) (a, T) + C(3) - J dyyln [l - ri^e^] J 



47ra 2 E 
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k B T 
lQna 2 

" dy 1 r dyy 2 
ev-l ^4 ev-l 



(35) 



where Fg (a, T) is the free energy in the plasma model given by Eq. (|24|), and r±(y) is defined in Eq. (|l9|). It is 
notable that the results of numerical calculations by this formula and by Eqs. (16), ( |30| ) (with different approaches 
to the zero-frequency term) coincide with an accuracy of 0.06% at a = 0.4 /xm and better than 0.01% for a > 3^m. 
Note, as discussed above, the Drude model leads to satisfactory s(i£) only up to £ ~ 10 15 rad/s and is in strong 
disagreement with the optical tabulated data for higher frequencies. Because of this, at T = 300 K it is meaningless 
to use the Drude dielectric function at separations < 0.4 /im. However, even at these separations, Eq. ( |35| ) is correct 
with an accuracy of 0.07%. 

As evident from Eqs. (|32|)-(|3"5|), the contribution of the zero-frequency term into the difference of the free energies 



AFp 



F D 

E 



Fg , computed in the framework of the Drude and plasma models, depends on the approach used: 
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(36) 
(37) 
(38) 



In the case of approaches a) and c) , the difference of the free energies contains linearly decreasing with temperature 
terms [Eqs. (^6|), (§)]• In the case of approach b), owing to the relaxation parameter, AFg falls off more quickly 



with decreasing temperature [the same is true for F^ , as is seen from Eq. (25), and for the summation term in the 
right-hand side of Eq. fl35|)]. It should be particularly emphasized that the presence of the linear terms in temperature 
in the free energy is in contradiction with the requirements of thermodynamics (see Sec. V). 

To obtain the approximate perturbative expression for the energy by analogy with Eq. (|3^), one should use the 
explicit dependence of 7 on temperature. It has been known that at temperature T > Td/4, where To is the Debye 
temperature (for Al T D = 428 K @), dj/dT = 7/T, i.e. 7 is linear in temperature. Generally, dj/dT = v-f/T with 
v = v{T) > 1. In Fig. 3, the dependence of 7 on temperature is plotted for Al on the basis of tabulated data p7[ . 
Finally, the required expression for the energy is 



E D (a,T) = E pl (a,T) + e { Q a ^ c) (a, T) + 



k B T 7 
47ra 2 lo v 



E 
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26 2 



e& - 1 
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dy 



e« - 1 



v - 1 



dyy 2 
ev - 1 



(39) 



where e$ = = 0, and 
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e °^ T) = -^- p (40) 

In Fig. 4 the results of the numerical calculations are shown for Al described by the Drude model in different 
approaches at T = 300 K. In the vertical axis the dimensionless ratios are plotted 

D E D (a 1 T) Fg(a,T) E D (a,0) 

\E D (a,Q)\' \E D (a,0)\' \E D {a,0)\ 1 ' 

as a function of the surface separation. Curve 1 shows the behavior of the relative energy (which is practically the 
same in all three approaches); curves 2a, 2b and 2c show the relative free energy in the approaches a), b), and c), 
respectively. The dashed curve is for the energy at zero temperature. All calculations are performed both using the 
exact expressions (|l|), ([H]) and the approximate ones (^), (S^) with coinciding results. 

It is important to explain in more detail the notation £/(a, 0). It is the value of energy in the framework of 
the Drude model (|30|), computed at zero temperature in the sense that Eq. ( pl| ) with a double integral instead of a 
discrete sum is employed. At the same time, in calculations of E D (a,0) the value of the relaxation parameter 7 at 
T = 300 K is used. We divide the calculational results into this hybrid quantity, previously used in literature (see, 
e.g., Refs. p|, p2] , p3| , p4|j2^ ] ) . This allows one to associate this quantity with energy and free energy in order to clarify 
its physical meaning. 

As is seen from Fig. 4, curve 1, illustrating the equal behavior of energy in all three approaches, and curve 2b, 
illustrating the behavior of free energy as given by the approach b), demonstrate plausible properties. Among other 
things, the free energy approaches energy with a decrease of the surface separation distance (compare with Fig. 2 in 
the case of the plasma model). As to the curves 2a and 2c, representing the free energy in the approaches a) and c), 
they do not approach to each other nor to energy within the application range of the Drude dielectric function. Note 
that even at separations of about 0.4 — 0.5 /im the free energy Fj? (curve 2a), obtained by the direct application of 
the Lifshitz formula, differs by 8% from the double integral E D (a,0) (dashed line). 

An important point is that not only the free energy of curve 2a but also 2b and 2c, and energy of curve 1 do not 
approach the dashed line in Fig. 4 representing the quantity which is in common use as a measure of energy at zero 
temperature §J||,||,||,||. This is clearly seen from Fig. 5 where the curves 1, 2b and 2c are reproduced on an 
enlarged scale for the smallest separations where the Drude model is applicable. The long-dashed curve 3 in Fig. 5 
illustrates the dependence of one more quantity on surface separation defined as 

n E°(a,T) 

where E® is the energy at a temperature T computed on the assumption that 7 does not depend on temperature 
(and preserves its value as at T = 300 K). Curve 3 is computed by Eq. (|3l|) with d"f/dT — 0. The same curve is 



obtained by the application of the approximate Eq. (39) with v = 0. 

From Fig. 5 we notice that curve 3 approaches the short-dashed curve with a decrease of a separation distance. 
Because of this, it may be concluded that the hybrid quantity E D (a,0) computed in the literature is in fact some 
approximation for E® , i.e. for the energy at temperature T computed without regard for the explicit dependence 
of the dielectric properties on temperature [remind that this kind dependence is absent in the case of the plasma 
model (see Sees. II, III) but is essential for metals described by the Drude model]. From Fig. 5 it follows that at a 
separation of 0.5 /im E D (a, 0) departs from the correct value of energy (curve 1) by approximately 0.75%. As to the 
free energies of the approaches a) and c), the deviations are larger (8% and 3.3%, respectively; these approaches are 
in contradiction with thermodynamics, see Sec. V). The above deviations should be added to the errors of E D (a,Q), 
discussed in Ref . [E3| , that are connected with uncertainties in the optical tabulated data. 



V. ENTROPY FOR THE THERMAL CASIMIR FORCE BETWEEN REAL METALS 

Considerations of the entropy of the fluctuating field in dependence on temperature allows one to test different 
approaches discussed above for conformity to thermodynamics. Entropy of the fluctuating electromagnetic field can 
be expressed in terms of a free energy 

dF E (a,T) 

S{a,T) = — (43) 
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or, taking into account Eq. (||), identically, as 



S(a,T) 



T 



[E(a, T) — Fe (a, T)} . 



(44) 



So it can be simply computed by the use of the results for the free energy and energy obtained in Sees. Ill, IV. 

Let us start with the plasma model where the analytical calculation is possible [approach d)] . At small separations 
(low temperatures) one can use Eqs. (p5[), ( |26| ) for the free energy and energy, respectively (a > X p is supposed). Then 
both Eqs. (43) and (H) lead to one and the same result 



S pl {a,T) = 



3fc B C(3) / T 
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(45) 



Note that this expression was first obtained in Ref. |34j with errors in numerical coefficients, because in Ref. gj] 
the energy of thermal photons was not taken properly into account. At large separations (high temperatures) the 
asymptotic expressions (^?]), (^8|) are applicable leading to 



(we have omitted exponentially small terms in 2irT/T e ff) 
of entropy for plates made of ideal metal 
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Eqs. @ 
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(46) 

]) lead to the values 
(47) 



for T <C T e ff, T ^> T e ff, respectively. The results (47) coincide with those obtained for an ideal metal in Ref. 
Asymptotical behavior of the entropy for an ideal metal in a high temperature limit was obtained also in Ref. 
The result of is, however, two times smaller than in Eq. ( |47| ) and Ref. [Q due to an error contained not only in the 
entropy but also in the expression for the Casimir energy between two plates made of ideal metal at zero temperature 
as is used in Ref. p9| |. 

It is obvious that Eq. ( f45| ) leads to nonnegative values of entropy with S pl (a, 0) = as is demanded by the third 
law of thermodynamics (the Nernst heat theorem j5(| ) . 

We now direct our attention to the entropy in the framework of the Drude model. As before, numerical calculations 
can be performed by the exact formulas for the energy and free energy or by the approximate ones valid at a > X p 
with coinciding results. From Eqs. (|35|), (|43| ) one obtains 



S D (a, T) = (a,T) + S^ c) (a,T) + ^-I^ 



i=i 



' -{» + 2)6 



dy 



ev - 1 



v 

Jul 



dyy 2 
e» - 1 



(48) 



Here 5 p! (a,T) is the entropy in the framework of the plasma model computed by Eqs. (p4|), (^3|), and SQ & ' h ' c \a,T), 
defined by 



(a,b,c) 




(a,T) 



i) 



167ra 2 dT 



ft' h ' c) (a, T) + C(3) - / dy y In (l - r± (yje"") 



(49) 



describes the contribution of the zero-frequency term of the Lifshitz formula into entropy in different approaches. 
Using the same perturbation expansions as in Sec. IV, one obtains 



S^(a,T) = 
S< h \a,T) = 
S ( o C) (a,T) = 



fcsC(3) 



2 \ 



XI 
3- p 



48a : 



(50) 



fe^C(3) X p 
87ra 2 7ra 



1 



3 X p 
2 7ra 



The results of numerical calculations using Eqs. (48), ( |50| ) for a = 2 /im are presented in Fig. 6. As is seen from the 
figure, in the approach a) entropy is negative in a wide temperature range from T = to almost T = 300 K, which is 
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a nonphysical result. In the approach a) entropy preserves the negative sign for lesser separations between the plates 
as well. In the approaches b), c) entropy is positive as it must be. In the approach b) S D (a,0) = 0, whereas in the 
approaches a), c) S D (a, 0)^0 which is in contradiction with the Nernst heat theorem. From Eqs. (ft8|), ( |50| ) it follows 
that 

S D (a,0) = S^ h ' c) (a,0), (51) 
^ b) (a,0) =0, S^(a,0) -4 a) (a,0) = 

where the absolute values of S , Q a '°' l (a, 0) are given by Eq. (|5(]). They are not only different from zero but depend 
on the parameters of the system (plate separation distance and plasma wavelength) which is prohibited by the third 
law of thermodynamics |30|. Because of this, approaches a) and c) must be rejected. Note also that approach a) 
predicts nonzero value of entropy at zero temperature for an ideal metal in contradiction with the field-theoretical 
result of Ref . 0] . As for approaches b) and d) , based on the special modification of the zero-frequency term of the 
Lifshitz formula and on the use of the plasma model, respectively, they are in agreement with the requirements of 
thermodynamics. To decide between them some additional considerations, which are presented in the next section, 
are needed. 



VI. CONCLUSIONS AND DISCUSSION 



In the above the correlation between the Casimir energy and the free energy at a temperature T is investigated 
for the case of two plane parallel plates made of real metal. It is shown that for the nondissipative media described 
by the real dielectric permittivity with no explicit dependence on temperature the photons between plates can be 
considered as the elementary excitations of the electromagnetic field interacting with a matter of plates. In this case 
the energy at temperature T is proved to be a sum of the (renormalized) energy of zero-point oscillations and thermal 
photons. If the media are dissipative and their dielectric permittivity depends on temperature, the simple picture 
above is not correct. The concept of thermal photons loses immediate significance and the energy of fluctuating field 
contains additional terms depending on the derivatives of the dielectric permittivity with respect to temperature. 

The expression for the energy at a temperature T found in this paper helps to elucidate the meaning of the so called 
"Casimir energy at zero temperature" calculated by many authors as a double integral using the Drude model and 
optical tabulated data at room temperature (note that this quantity is of great importance as it is proportional to the 
Casimir force in the configuration of a sphere or a spherical lens above a plate used in experiments [pp^flM^lUiaJla] ) . 
The commonly accepted opinion that the above-mentioned quantity is approximately equal to the free energy at small 
temperatures (small separations) is inexact. In fact, even at rather small separations the "Casimir energy at zero 
temperature" deviates from the free energy by several percent but approaches to the energy at room temperature 
calculated on the assumption that the dielectric permittivity does not depend on T explicitly (this assumption is not 
correct in the case of the Drude dielectric function) . 

Different approaches to describe the thermal Casimir force from recent literature were compared and analysed 
[approaches a), b), c) in the framework of the Drude model and approach d) in the framework of the plasma model - 
see Introduction]. The quantitative expressions for the entropy of the fluctuating field are obtained here for the first 
time in the case of real metals. They give the possibility to conclude that the approaches a) and c) are in contradiction 
with the principles of thermodynamics and must be rejected. The approaches b) and d) are found to be in agreement 
with thermodynamics. 

To make a choice between the approaches b) and d) let us discuss the behavior of the dielectric permittivities 
of the plasma and Drude models at small frequencies. Several authors [^ , [30|]3l|]35fl give preference to the Drude 
model because it shows uj^ 1 frequency dependence of the dielectric permittivity at small frequencies as it follows from 
Maxwell equations (compare with u>~ 2 frequency dependence given by the plasma model). Although this statement 
is true, it should be remembered that the Drude model is not applicable at all frequencies. We note that the concept 
of e(oj) itself, not only the Drude model, does not work in the domain of the anomalous skin effect (see Sec. IV). 
As to the quasistatic limit, although s(u>) is of order w _1 in this domain, the Drude model is also not applicable as 
the correct e(to) is pure imaginary. Since the zero-frequency term of the Lifshitz formula necessarily belongs to the 
domain of the quasistatic fields, where the concept of traveling waves fails, the substitution of the Drude dielectric 
function into this term (resulting in all the above problems) seems to be unjustified. 

To clarify the situation with the thermal Casimir force, let us consider two pairs of plane parallel plates a — 5 /im 
apart made of Al (one pair of plates) and of indium tin oxide (the other one). Due to large separation distance, the 
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asymptotic of high temperatures is applicable and only the zero-frequency term of the Lifshitz formula determines 
the total value of the Casimir force. At quasistatic frequencies both Al and indium tin oxide are good conductors. 
Because of this, the Lifshitz formula would lead to one and the same Casimir force at 5 /xm separation for both pairs 
of plates if one substitutes into it the actual reflection properties of these materials at zero frequency. This conclusion 
is in contradiction with intuition. Note that an indium tin oxide is transparent to visible and near infrared light. 
Within a wide wavelength range 7/xm< A < 100 /mi around the characteristic wavelength A c = 62.8 /jm (the latter 
corresponds to the characteristic frequency lu c — c/2a = 3x 10 13 rad/s), giving the main contribution into the Casimir 
force at zero temperature || , the reflectivity of indium tin oxide is below 80% |5l| . Note that the second parameter 
of the problem, first Matsubara frequency, is % = 2itkBT/Ti = 2.45 x 10 14 rad/s, i.e. Aj\/ = 7.7 /tm, which belongs to 
the region of even larger transparency of indium tin oxide. In this situation it is difficult to imagine that at a = 5 /iin 
the indium tin oxide plates are attracted with the same Casimir force as Al plates which are almost perfect reflectors 
within a wide range around the characteristic wavelength. 

We can avoid this contradiction between the literally understood theory and physical intuition if we assume that it is 
not correct to substitute the actual behavior of the dielectric permittivity at zero frequency into the Lifshitz formula. 
Instead, in order to obtain the physically correct results, the frequency dependence of the dielectric permittivity 
and reflection coefficients around the characteristic frequency should be extrapolated to zero Matsubara frequency 
and substituted into the Lifshitz formula. If this conjecture is accepted, one should conclude that within the range 
of micrometer separation distances between plates the plasma model dielectric function, i.e. the approach d), is 
preferable as compared with the use of the Drude dielectric function combined with any of the above approaches a), 
b), c). It is apparent from the fact that the plasma dielectric function and respective reflection coefficients admit 
reasonable continuation from the range of infrared optics to zero frequency. 

The contradictions discussed in this paper lead to a conclusion that the concepts of the frequency dependent 
dielectric permittivity and fluctuating electromagnetic field inside media in application to the thermal Casimir force 
between real metals are inadequate idealizations. Less sophisticated approaches, such as the surface impedance 
approach (the Leontovich boundary conditions), which does not consider the fluctuating field inside matter [j3|,j33|, 
appear to be more adequate and lead to physically justified results for all separation distances between plates. By 
way of example, in the domain of the infrared optics the surface impedance leads to the same results as the Lifshitz 
formula in combination with the plasma model [approach d)] . If the characteristic frequency belongs to the domain 
of the normal skin effect, where the Drude model is physically correct, there is no reasonable continuation of e to zero 
frequency avoiding the above problems. At the same time, the impedance approach, when applied in the domain of 
the normal skin effect, leads to quite satisfactory results j3^[ coinciding with those for ideal metal as it must be at 
separations larger than 0.1cm [almost the same results are given in this domain by the approaches b), c)]. 

To conclude, at present the impedance approach can be considered as the most universal, reliable and straightforward 
way to calculate the thermal Casimir force between real metals at different separation distances. In the domain of 
micrometer separations the plasma model is also realistic. Regarding the Drude model, it can be used to describe 
the thermal Casimir force only with some appropriate modification of the zero-frequency term of the Lifshitz formula 
[like in the approach b), for instance]. 
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FIG. 2. Relative energy at temperature T = 300 K (curve 1), free energy (curve 2), and energy at zero temperature (dashed 
line) versus surface separation in the framework of the plasma model. 
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FIG. 4. Relative energy at temperature T = 300 K (curve 1), free energy [curve 2a in the approach a), curve 2b in the 
approach b), and curve 2c in the approach c)], and energy at zero temperature (dashed line) versus surface separation in the 
framework of the Drude model. 
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FIG. 5. Relative energy at temperature T = 300 K (curve 1), free energy [curve 2b in the approach b) and curve 2c in the 
approach c)], and energy at zero temperature (short-dashed line) versus surface separation in the framework of the Drude model 
reproduced on an enlarged scale. Long-dashed curve 3 presents energy computed on the assumption that dielectric permittivity 
does not depend on temperature. 
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FIG. 6. Entropy of fluctuating electromagnetic field in the framework of the Drude model versus temperature computed on 
the basis of approaches a), b), and c) (curves a, b and c, respectively). 
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